A simple variational approach to the 
quantum Frenkel-Kontorova model 

Choon-Lin Ho^'^ and Chung-I Chou^ 
^Department of Physics, Tamkang University, Tamsui 25137, Taiwan* 
'^Theory Division, Institute of Particles and Nuclear Studies, KEK, 
Tsukuba, Ibaraki 305, Japan 
^Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 

Abstract 

We present a simple and complete variational approach to the one- 
dimensional quantum Frenkel-Kontorova model. Dirac's time-dependent vari- 
ational principle is adopted together with a Hatree-type many-body trial 
wavefunction for the atoms. The single-particle state is assumed to have 
the Jackiw-Kerman form. We obtain an effective classical Hamiltonian for 
the system which is simple enough for a complete numerical solution for the 
static ground state of the model. Numerical results show that our simple ap- 
proach captures the essence of the quantum effects first observed in quantum 
Monte Carlo studies. 
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1. Introduction. The Frenkel-Kontorova (FK) model |T],0 is a simple one-dimensional 
model used to study incommensurate structures appearing in many condensed-matter sys- 
tems, such as charge-density waves, magnetic spirals, and adsorbed monolayers 0|. These 
modulated structures arise as a result of the competition between two or more length scales. 
The FK model describes a chain of atoms connected by harmonic springs subjected to an 
external sinusoidal potential. In an important development in the study of the classical 
FK model, Aubry @] first made use of the connection between the FK model, the so-called 
"standard map" , and the Kolmogorov- Arnold- Moser (KAM) theorem to reveal many inter- 
esting features of the FK model. Particularly, he showed that when the mean distance (also 
called the winding number) between two successive atoms is rational, the system is always 
pinned. But when the winding number is irrational, there exits a critical external field 
strength below (above) which the system is unpinned (pinned). This transition is called by 
Aubry a "transition by breaking of analyticity" , and is closely connected with the breakup 
of a KAM torus. It is very analogous to a phase transition, and various critical exponents 
and questions of universality have been extensively studied in the past. 

In recent years, the FK model has been applied to the study of transmission in Josephson 
junction and atomic-scale friction-nanoscale tribology In these cases, quantum effects 
are very important. Unlike the classical case, study of quantum FK models is rather scanty. 
It was first considered in a quantum Monte Carlo (QMC) analysis in 0]. Their main obser- 
vation is that the map appropriate to describe the quantum case is no longer the standard 
map, but rather a map with a sawtooth shape. An explanation of this phenomenon was 
later attempted in using a mean field theory in which the inclusion of the contribution 
from quasidegenerate states is essential. But the mean field computations with these states 
are rather involved, and the quantum map that is to replace the classical standard map in 
different quantum regimes is not clearly identified. More recently, a less complicated ap- 
proach was proposed in which uses the squeezed state function to demonstrate that the 
sawtooth behavior is simply the result of quantum fluctuations. In our opinion, the approach 
adopted in P] is very appealing in principle. However, we believe that some difficulties in 
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this work need to be overcome before it could be considered satisfactory. First, the assumed 
squeezed state many-body ground state is general enough so as to include the correlations 
of the positions of the atoms, expressed by the covariances Gij = {{xi — Xi){xj — Xj)) {i ^ j), 
where Xi is the position of the ith atom, (■ ■ ■) is the expectation value in a given quantum 
state, and = (xj). However, to find the equilibrium state of the model, one has to solve 
a system of coupled equations of the variables Xi and the Gij. The equations obtained are 
so complicated that the task of solving them within a single numerical framework is very 
difficult. In fact, in a hybrid numerical analysis was adopted in which the equations 
for the Gij were not solved. Instead, the values of Gij were taken from QMC data. These 
values were then treated as initial conditions in solving the equations for the atomic posi- 
tions Xj's. Technically, such hybrid analysis is not satisfactory. Second, the covariance terms 
Gij {i 7^ j) are constrained by the values of the fluctuation terms Ga and Gjj through the 
Cauchy-Schwarz inequality. These constraints guarantee the boundedness from below of the 
effective Hamiltonian ||^. But then this also calls for a proper variational principle that has 
to take care of the interdependence of the Gij terms. 

In this letter we shall show that all the essential features observed in the QMC studies can 
be obtained from an independent-particle picture of the many-body ground state without 
the covariance terms. In the independent-particle picture the many-body trial wavefunction 
are factorizable into single-particle states. One can assume the single-particle state to have 
the form of a squeeze state. For the quantum FK model, a simpler and, in our view, more 
elegant approach is to use the Jackiw-Kerman (JK) function as the single particle state. 



We shall show that this simple independent-particle approach produces an effective classical 
Hamiltonian which is bounded below, admits simple numerical solution of the ground state 
without recourse to QMC analysis, and reproduces the essential features observed in QMC 
studies. 

2. Effective Hamiltonian. The Hamiltonian of the quantum FK model is given by 
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Here qi and pi are the position and momentum operators, respectively, of the zth atom, 7 the 
elastic constant of the spring, V and 2t[/Iq are the strength and the period of the external 
potential. As in ||^, it is convenient to use the dimensionless variables Qi = loqi, Pi = 
^oPi/ a/?ti7, and K = Vll/j. With these new variables, we obtain the following dimensionless 
Hamiltonian H 
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We have H = 'jH/Iq. The effective Planck constant is h = Ml/ ,ynvj. For the classical FK 
model, the Aubry transition occurs at the critical value Kc = 0.971635 ■ • -. 

To study the ground state properties of the quantum FK model in (^, we adopt here 
the time-dependent variational principle pioneered by Dirac |Tl|. In this approach, one first 
constructs the effective action T = J dt {^,t\iTidt — 1-L\^,t) for a given system described by 
Ti and t). Variation of F is then the quantum analogue of the Hamilton's principle. The 
time-dependent Hatree-Fock approximation emerges when a specific ansatz is made for the 
state \^it). We now assume the trial wavefunction of the ground state of our quantum FK 
system to have the Hatree form t) = Hi iV'i; ^)) where the normalized single-particle state 
is taken to be the JK wavefunction |T0| : 
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The real quantities Pi{t), Gi{t) and Hj(t) are variational parameters the variations 

of which at t = ±00 are assumed to vanish. The JK wavefunction can be viewed as the 



Q-representation of the squeeze state |jT2|. We prefer to use the JK form since the physical 
meanings of the variational parameters contained in the JK wavefunction are most trans- 
parent, as we shall show below. Furthermore, the JK form is in the general Gaussian form 
so that integrations are most easily performed. 

It is not hard to check that Xi and pi are the expectation values of the operators Qi and 
A: Xi = Pi = (^lAl^). Also, one has - = hG^, and {^\thdt\^) = 

HiiPi^i — ^Gjllj), where the dot represents derivative with respect to (w.r.t.) time t. It is 



now clear that hGi is the mean fluctuation of the position of the i-th atom, and that Gi > 0. 
With these expectation values, the (rescaled) effective action F for the dimensionless H can 



be worked out to be T{x,p, G, 11) = / dt [J2i^o^ (Piii + ^HiGj) — Hejf], where luq = y'y/m 
is the angular frequency of the spring, and i^e// = (^l-^l^) is the effective Hamiltonian 
given by 

+ E ^ (■^'•+1 -^if 
i ^ 

+ E I (G'.+i + G.) 
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- ^i^exp cosxj . (4) 

The last term in 

(I) can be very easily obtained from {%F{Qi)\^) = E^^o F^'^"'\xi){hGi)"' /{2m)U, where 
F^"'\x) = d"'F(x)/dx^, and n\\ = n{n — 2)(n — A) • • ■ 1. Eq. (^) is bounded from below. One 
sees from the form of the effective action F that Flj is the canonical conjugate of Gj. 

Varying F w.r.t. x, p, G and FI then gives the equations of motion in the Hatree-Fock 
approximation. Since we are mainly concerned with the static properties of the ground state 
of the quantum FK model, we must set the time derivatives of these variables to zero. This 
gives the equations which determine the values of variational parameters corresponding to 
the equilibrium states (which include the ground state). Equivalently, we can obtain the 
equations for the equilibrium states by directly varying the effective Hamiltonian Heff w.r.t. 
the variables. Varying Heff w.r.t. pi, Flj, Xi and Gi give, respectively, 

Pi = 0, 4n,Gj = , (5) 
x.,.-2.. + .._,^^irexp(-^G.)sm.., (6) 

^Gr2 _ Kexp (^'^'^i^ cosx^ - 2 = AUj . (7) 

The second equation in implies Flj = as > 0. This in turn means that the right 
hand side of eq.P) is equal to zero: 



- K exp (^-^'^^ cos(xi) -2 = 0. (8) 

In the limit h = 0, eq.(|^) is equivalent to the standard map. We note that eq.(|^) was also 
obtained in §]. This is because in the formulation in the covariances Gi/s decoupled from 
the Xi and the fluctuations Gu {Gi in our case) in the variation of their Hamiltonian w.r.t 
Xj. Unlike our case, of course, these covariance terms do actually influence the solutions of 
(^) through other equations obtained by variation of the Hamiltonian w.r.t. the Gu and 
Gij. And it is these equations that caused the difficulties mentioned in the Introduction. In 
particular, the values of the Gu were input from the QMC data in order to solve for the Xi 
in (P). Our simple approach, on the other hand, allows us to solve for both the values of Xi 
and Gi coupled by eqs.(||) and (H) consistently by a single numerical method. 

^From pj = Ilj = and (H), we see that the problem of finding the static ground state 
of the quantum FK model reduces to the problem of minimizing w.r.t. to Xi and Gi the 
following effective potential 
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Eq.(^ and (|) are just the conditions dVejf/dxi = and dVefj/dGi = 0, respectively. 

3. Numerical results. We numerically solve for the set of variables Xi and Gi which char- 
acterize the ground state using the Newton method In all our numerical computations 
the winding number P/Q = 610/987, which is an approximation of the golden mean winding 
number (a/S — l)/2, is used with the periodic boundary condition Xi+Q = Xi + 2txP. This 
winding number is much more accurate than those used in previous works to approximate 
the golden mean number, thus giving us better accuracy in the computations of physical 
quantities related to the ground state. We emphasize that all values of Xi and Gi are deter- 
mined by the same numerical method consistently. In particular, we do not have to input 



the values of Gi from quantum Monte Carlo results in order to solve for the Xi. 

Having obtained the values of Xi which give the mean positions of the quantum atoms in 
the chain, we can compare the results with the classical configuration, following p|, in two 
ways: (1) by the quantum hull function, which is the plot of Xi (mod 2tt) of the atoms against 
their unperturbed positions 2niP/Q (mod 2n); (2) by the so-called ^f-function, defined by 

gi = K'^ (xi+i - 2xi + Xi-i) (10) 

versus the actual atomic positions Xj. /^From (|^), we also have 

9i = exp (^-^G^j sinxi . (11) 

Here Gi is related to Xi by eq.(||). We see from this equation that quantum fluctuations Gi 
will modify the shape of the classical sine-map. In addition to these two types of graphs, 
we also plot the graph of Gi against the unperturbed and the actual positions. The formal 
graph was first introduced in to show the strong correlation of the fluctuations of atoms' 
positions with their unperturbed positions. We introduce the latter type of graphs here 
since we think that such graphs provide better picture about how the quantum fluctuations 
of the atoms are related to their actual positions. 

In Fig. 1 we show the four graphs mentioned above with different values of h for the 
supercritical case K = 5. Fig. 1(a) shows the quantum hull functions. For small values of 
h the quantum hull function consists of a countable set of steps discontinuities, just as in 
the classical case: the atoms are in a pinning phase. In fact, the atoms are more likely to 
be located near the valley of the external potential well, namely, near Xi = (mod 2it). 
As the quantum effect increases, i.e., for increasing values of h, the quantum hull function 
gradually changes into a monotonic analytic function, signifying that the system is entering 
the depinning phase. There exists a critical value, approximately = 6.58 for K = 5, above 
which the quantum hull function changes from an nonanalytic function to an analytic one. 
This is a quantum analogue of the Aubry transition in the classical case, and can therefore 
be called the quantum Aubry transition. 
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Next in Fig. 1(b) we show the graphs of the (yf-function. The curve defined by ( [TlD 
with Gi satisfying are shown here as dashed curves for different h. In the classical limit 
{h = 0) this curve is simply the standard map (sine-curve). As h increases, the amplitude of 
the curve decreases. For sufficiently large h, the curve resembles more closely a "sawtooth" 
shape. This is first noted in QMC study in . Here we see that it comes out very naturally 
from the equation of motion (||) and (|ll]). In the supercritical case {K = 5), when h < he, 
the positions Xi of the atoms cover only a subset of the (^-curves. This is in accord with the 
fact that the atoms are in the pinning phase [cf. Fig. 1(a)]. As h increases, the points begin 
to spread along the g-cmve. When h > he, the (7-graph is completely covered as the system 
has entered the depinning phase. 

Fig. 1(c) shows the quantum fiuctuations Gi plotted against the actual atomic positions 
Xi. The dashed curves represent the curves of eq.(P) for different h. For small h, the atoms 
are located near Xi = (mod 2tt) with small values of Gi which means, from (Q), that 
the wavefunctions are highly peaked at these positions. As the quantum effect increases, 
the external potential is so modified that now the atoms could be found at other positions, 
but with atoms at Xj = vr (mod 2tt) having the largest value of Gj. This indicates that 
wavefunctions of the atoms near the top of the potential are more extended with smaller 
amplitudes. Again, when h > he, the curves of (|^) are completely covered by the solutions Xj. 
To compare with the results in , we plot the values of Gi against the unperturbed positions 
in Fig. 1(d). One sees that the values of Gi are strongly correlated with the unperturbed 
positions, as first noted in 0. For h < he the graphs consists of steps discontinuities, and 
for h > he the graphs are continuous. This is correlated with the graphs of the quantum 
hull function in Fig. 1(a), since from (H) any fixed value of Xi correspond to a fixed value of 

Next we show in Fig. 2 the corresponding graphs for the case K = 1.5. This represents the 
situation which is slightly over the critical classical case. The general trends of the behavior 
of the graphs are the same as those in Fig. 1. As expected, quantum Aubry transition takes 
place at a smaller he = 1.17. We note here that the shape of the ^f-function at large h in 

8 



this case is intermediate between a sine and a sawtooth map. 

We have also checked the subcritical cases with K < K^.. The classical system is already 
in the depinning phase in this regime. Quantum fluctuations only enhance the trend of 
depinning. The gf-function is found to be closer to a sine-shape with smaller amplitude for 
higher h. This is consistent with the QMC results Q. 

Finally, we note here that, while we have reproduced the essential features first observed 
in the QMC studies of the quantum FK model, there is also slight discrepancy between the 
results of these two approaches. The difference is that, for a fixed value of the QMC 
results indicated that the sawtooth shape of the (7-function appeared at a lower value of 

and that the atoms began to spread along the (yf-curve also at a smaller h. For example, 
at = 5 the QMC results showed that the above situation already appeared at ^ = 0.2, 
while our results (cf. Fig. 1(b)) indicate that at = 5 and at a higher h = 2 the system 
is still closer to the classical case. We believe this could be explained as follows. First, our 
independent-particle wavefunction is only the lowest order approximation of the many-body 
wavefunction of the quantum FK system. A more accurate description of the system will 
require a better assumption of the wavefunction than that assumed here. This presumably 
may require the inclusion of the effects of the covariance terms as advocated in |Q, but with a 
more appropriate variational principle to circumvent the difficulties already mentioned in the 
Introduction. Second, our results are obtained at zero temperature, while those in the QMC 
analysis were obtained, by the nature of the method itself, at small but finite temperatures 
(temperature T=0.0067 as given in 0). It is natural that thermal fluctuations will cause 
the atoms to spread away from their zero-temperature positions. 

4- Summary. In conclusion, we have presented a simple and complete variational ap- 
proach to the quantum FK model based on a Hatree-type many-body trial wavefunction of 
the JK form. The effective Hamiltonian obtained is bounded below, and is simple enough 
for a complete numerical solution for the static ground state of the model in various quan- 
tum regimes. Numerical results show that our simple approach captures the essence of the 
quantum effects first observed in QMC studies. The map appropriate for the quantum FK 
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model is well described by eq.(|5) and (||). In contrast to previous approaches, we do not 
require the existence of the comphcated quasidegenerate states, or the partial help from 
quantum Monte Carlo data in order to obtain these results. 
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Figures captions 

Fig. 1 Structure of the quantum ground state for K = 5 and winding number P/Q = 
610/987 at h = 2 (black dots), 6 (white dots) and 7 (black curve), (a) quantum hull 
function plotted against unperturbed atomic positions; (b) (^-function plotted against 
actual atomic positions (the dashed curves represent eg. (pA]) with Gi satisfying (H); (c) 
and (d) quantum fluctuations Gi plotted against the actual and unperturbed positions, 
respectively. The dashed curves in (c) represent the curves of eq.(P) for different h. 

Fig. 2 Same as Fig. 1 for i^T = 1.5 and h = 0.5 (black dots), 1.0 (white dots) and 2 (black 
curve) . 
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